Sound propagation in elongated superfluid fermion clouds 
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We use hydro-dynamic equations to study sound propagation in a superfluid Fermi gas inside a 
strongly elongated cigar-shaped trap, with main attention to the transition from the BCS to the 
unitary regime. We treat first the role of the radial density profile in the quasi-onedimensional limit 
and then evaluate numerically the effect of the axial confinement in a configuration in which a hole 
is present in the gas density at the center of the trap. We find that in a strongly elongated trap 
the speed of sound in both the BCS and the unitary regime differs by a factor \/3/5 from that 
in a homogeneous three-dimensional superfluid. The predictions of the theory could be tested by 
measurements of sound-wave propagation in a set-up such as that exploited by M.R. Andrews et al. 
[Phys. Rev. Lett. 79, 553 (1997)] for an atomic Bose-Einstein condensate. 
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Strong evidence for a superfluid state in ultracold 
Fermi gases near a Feshbach resonance has come from 
experimental studies using as probes ballistic expansion 
i| [J 0] collective modes 0, Q , RF spectroscopy 
3, and the generation of quantized vortices Since 
the BCS weak-coupling limit and the unitary strong- 
coupling limit are characterized by the same collective- 
mode frequencies and expansion dynamics 0, atten- 
tion has been drawn to the study of sound-wave prop- 
agation as allowing a clear identification of these two 
regimes. In particular Ho 0] and Heiselberg have 
evaluated the speed of first and second sound in a ho- 
mogeneous Fermi superfluid as functions of the coupling 
regime. At low temperature the first sound velocity is 
given by U\ — + (2/TT)hkp a)/3 ~ in the 

dilute BCS limit and by u\ ~ 0.37?;f in the unitary limit. 
Here vp and k-p are the Fermi velocity and wavenumber, 
while a is the s-wave scattering length. 

The main purpose of this Letter is to evaluate the prop- 
agation of density perturbations in a superfluid Fermi gas 
as a function of coupling strength in an experimentally 
realizable set-up. This simulates the set-up used in the 
experiments of Andrews et al. |13j on the propagation 
of sound pulses along the axis of an elongated cloud of 
Bose-Einstein condensed 23 Na atoms. In their experi- 
ments a density perturbation is generated by turning on 
a laser at the center of the trap and two pulses propagate 
in opposite directions along the trap axis. We describe 
the dynamics of density fluctuations in the superfluid by 
hydrodynamic equations, namely the continuity equation 



dtn + V • (nv) = 
and the Euler equation 
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for the time-dependent density profile n(r, t) and veloc- 
ity field v(r, t) of the trapped gas. The equation of 
state enters through the density-dependent chemical po- 
tential fj,(n), and V(r) describes the external potentials. 



The laser beam is simulated in our study by a time- 
independent effective potential U(z) — Uoe~ z ' w hav- 
ing amplitude Uq and width w, which is turned on at the 
center of the trap. 

The Q1D configuration. We discuss first the quasi- 
onedimensional (Q1D) case corresponding to V(r) = 
mw i- r i_/2i where u>± is the radial trap frequency. This 
is amenable to analytical treatment. The equilibrium 
density n<j(r±) is determined from the relation /x(no) = 
p, — mu)j_rj_/2, where \x is determined by the nor- 
malization condition J no(r±)d 2 r± = h, n being the 
linear density. Linearization around equilibrium sets 
n(r,i) = n + Sn(r ± )e i( - ut - qz ^ and fi(n) = n(n ) + 

dv/dn\ n =n§ n ( r t) el{uJt ~ qz) ■ Thus E( l s - and © lead 
to the eigenvalue equation 

muj 2 8n — q 2 (no dfi/dn\ n=n 5n) 

-V±-[n Q Vx(d t i/dn\ n= Jn)] (3) 

and integration in the (x, y) plane yields the dispersion 
relation 



= q (^—J n d(i/dn\ n = no 6nd 2 r±/ J Snd 2 r±^j 
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The lowest collective mode of the gas is obtained by tak- 
ing Sn as the lowest-energy solution of Eq. © with q = 0, 

i.e. 5n(r±) cx (d/i/dn\ n=ng ) . Therefore the sound ve- 
locity is 



(1) t»i= (J^J n d 2 r ± / J (dn/dn\ n=no ) 1 d 2 r A ^j 
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Equation (JSJ applies to all Q1D gases which can be de- 
scribed by hydrodynamic equations. 

For a power-law form of the equation of state, n(n) = 
Cn 1 , the equilibrium density profile takes the explicit 

form n (r±) = {fi/C) lh (l - r^/R^) 1 ' 1 where R TF = 
y / 2p,/moj 2 L is the radial size of the cloud. Equation (0) 
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then reduces to 
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For 7 = 1 Eq. JHJ) gives the speed of sound pulses prop- 
agating along the axis of a Q1D condensate of bosonic 
atoms or °f molecular bosons formed by pairs of 
fermionic atoms. For a Fermi gas in the BCS (fcp|a| <C 1) 
and the unitary (/cfH 2> 1) regime the equation of state 
is instead proportional to the Fermi energy ^-(37r 2 n) 2 / 3 , 
so that 7 = 2/3 and u\ = yj^fl/^m. In the former limit 
at very low temperatures p, = mv F /2 and u\ = vf/V§, 
while in the latter /2 = m(l + (3)vp/2 with [3 being a 
constant that has been evaluated in various approaches 
15]. Quantum Monte Carlo calculations find (3 ~ —0.57 
lij|l7j and then u\ ~ 0.29uf- These velocities thus differ 
by a factor ^/3/5 from those in the homogeneous Fermi 
gas |l2j and would be measured in a sound propaga- 
tion experiment through the central portion of a strongly 
elongated trap. It is worth noticing that these results 
differ from what a ID calculation would predict: in the 
latter it is straightforward to show that u\ = ^J^fi/m, 
i.e. u\ = i>f for a strictly ID hydrodynamic Fermi gas. 

At intermediate couplings, when k-p\a\ is finite, C de- 
pends on the density and several forms of the equation of 

I 1 1 1 1 1 

state have been proposed |15|j Ha 113 • A possible solution 
is to use again a power law with an effective exponent 7, 
which may be inferred from a measured collective mode 
frequency 01 • The effective value of C can instead be 
determined from the cloud size. In this work we focus 
on a gas with a negative scattering length and choose 
the equation of state given by Heiselberg in the Galitskii 
approximation This is 



M(n) = ^(3- 2 «) 2/3 
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with ci = 2a/(97r) 1 / 3 giving the mean-field contribution 
and c 2 = 6[11 - 21n2](37r 2 ) 1 / 3 a/(357r) the second-order 
contribution from the ladder diagrams, which takes into 
account the opening of the superfluid gap. Equation Q 
reproduces the result fi = m[l + 4aA:F/(37r)] v F /2 in the 
limit of a dilute system, and is approximately valid at 
intermediate densities in the unitary limit where it gives 
m(l + f3)v F /2. The ratio between c\ and c 2 in this case 
yields (3 ~ —0.67. The results for the sound velocity as 
obtained from a self-consistent evaluation of the chemi- 
cal potential and of the equilibrium density are shown in 
Fig. n In this calculation we have taken a linear den- 
sity n corresponding to the chemical potential of 2 x 10 3 
fcrmions in a cigar-shaped confinement with trap fre- 
quencies u!± = 100 s -1 and uo z = 6.3s" 1 (see below). 
The BCS and unitary limits are recovered for a — > 0~ 
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FIG. 1: Speed of sound ui in a Q1D Fermi gas (in units of 
the Fermi velocity uf) as a function of the s-wave scattering 
length a (in units of the Bohr radius ao and on a log scale). 
The full and the dashed line are the results from the equation 
of state 10 and from the mean-field (MF) model, respectively. 
The dot-dashed lines mark the BCS (upper) and the unitary 
(lower) limits. 



and a — > — 00, respectively. The dashed line in Fig. Q] 
shows the sound velocity obtained within the mean-field 
model for a dilute Fermi gas by setting c 2 = in Eq. 
(J2J). The mean-field model well describes the BCS limit, 
but as the coupling increases it predicts a strong drop 
of the sound velocity heralding collapse. The location 
and width of the transition from the BCS to the unitary 
regime are determined by the system parameters. 

Effect of axial potentials. A weak harmonic confine- 
ment along the z direction induces a smooth inhomo- 
geneity in the axial density profile, which can be treated 
within a local density approximation. Assuming again a 
power-law equation of state, Eq. JSJl leads to the velocity 
field 



vi(z) = 



H 7 



1/2 
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with Z> 



TF 



m 1 + 7 

2 , which corresponds to an oscil- 
latory motion with frequency uj = uj z / \[2~- ~ 2/7 and ve- 
locity u\ at the center of the trap. Even though u\ de- 
pends on C, the oscillation frequency is the same in the 
BCS and in the unitary limit, since the product between 
the wavenumber q cx 1/Ztf °^ vVA an d t ne velocity 
Mi oc %f\± does not depend on the chemical potential. 

We next carry out a numerical study of sound propa- 
gation using the hydrodynamic equations and J5J in 
a set-up simulating that used in the experiments of An- 
drews et al. [I3j. In addition to radial and axial confine- 
ments provided by a strongly elongated harmonic trap, 
we switch on a hole-digging potential U(z) at the cen- 
ter of the trap. The atoms are thereby expelled into 
two density perturbations that travel towards the ends 
of the trap. Since the perturbation velocity is mainly in 
the z direction and imaging techniques will provide col- 
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FIG. 2: Integrated density profile n (in arbitrary units) as a 
function of axial coordinate z (in units of a± = {%/ (majx)) 1//2 ) 
at different times starting from the bottom at t = and 
increasingly separated by 1.5/ui±. The left, middle, and right 
panel correspond to a/ao = — 80, — 10 5 , and — 10 6 . For the 
sake of visibility the plotting range is smaller than the full 
extent of the density profile. 

umn densities, we focus on the evolution of the radially- 
integrated density profile h(z). The numerical calcula- 
tions are carried out for a system of N = 2 x 10 3 40 K 
atoms equally distributed in two spin-polarized states, 
with trap frequencies lo± = 100s _1 and lo z — 6.3 s -1 . 
The potential U(z) has a width w equal to 5% of the 
axial width of the density profile and an amplitude Uo 
of magnitude up to hcu z . In Fig. [21 we illustrate the time 
evolution of h(z) for three values of the scattering length, 
a = -80, -10 5 , and - 10 6 Bohr radii. The size of the 
cloud decreases on increasing the coupling strength and 
at the same time the speed of the density perturbation 
also decreases. 

From the evolution of the integrated density we cal- 
culate the mean velocity (v z ) of the perturbation in a 
time interval At = The results for (v z ) as a 

function of the coupling strength and for two ampli- 
tudes of the perturbing potential are shown in Fig. [31 
The comparison is also illustrated between the case with 
A = uj z /uj± ~ 0.06 and that with A = 0.01 at the same 
zero-coupling Fermi velocity. For weak perturbations (for 
example Uq = 0.167kj 2 , full symbols in Fig. [SI, giving Sn 
of the order of a few percent of the equilibrium density at 
the trap center, the speed of sound for both values of A is 
found to be in good agreement with the theoretical pre- 
diction given in Eq. JSJ for the Q1D gas (thick solid line). 
The slight difference seen in Fig. [3] may be accounted for 
by the effect of the weak confinement causing the speed of 
sound to decrease while the density perturbation moves 
away from the central region of the trap (see Eq. JSJ). 
Beyond the perturbative regime (empty symbols in Fig. 
[3J referring to Uq = 0.8hio z and Sn/n ~ 20%) we find 
an upward shift of the mean velocity as a result of the 
increase in the local density. The overall behavior of (v z ) 
is still the same as in the perturbative regime. The ratio 
between the values at zero and large \a\ is about 0.61, to 
be compared with y/1 + (3 ~ 0.57 for both the Q1D gas 
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FIG. 3: Sound velocity in a superfluid Fermi gas as a func- 
tion of \a\ (in units of the Bohr radius ao and on a log scale). 
Circles and squares correspond to the average speed of per- 
turbations propagating in a trap with A ~ 0.06 and A = 0.01. 
Empty and full symbols refer to Uo/hui z = 0.8 and 0.16, re- 
spectively. The thick solid line is from Eq. 0, and the dot- 
dashed lines mark the BCS (upper) and unitary (lower) limits. 
The thin solid lines are guides to the eye. 

(A = 0) and the homogeneous gas. 

Finally, we compare in Fig. 0] the sound velocity for 
A ~ 0.06 in a superfluid Fermi gas with that in a normal 
Fermi gas at a temperature T = 0.2 Tp, T-p being the 
Fermi temperature at zero coupling. The dynamics of 
the normal gas has been studied by the Vlasov-Landau 
kinetic equations as described in Ref. [2(J. Our method 
involves binary collisions between the particles and in- 
cludes Pauli-blocking effects. At low coupling the normal 
Fermi mixture is in the collisionless zero-sound regime 
and sound propagation is associated with an anisotropic 
deformation of the Fermi surface, the zero-sound veloc- 
ity being close to vp as shown in Fig. 0] With increasing 
attractions the velocity of the perturbation drops. This 
drop occurs at about a ~ — 10 4 ao, and can be attributed 
to both the transition towards the first-sound velocity 
Vf I y/5 and the instability of the approaching collapse of 
the normal state. The transition to the superfluid state 
restores stability. 

In summary, we have studied sound propagation in a 
superfluid Fermi gas confined inside a strongly elongated 
harmonic trap. We have found that the sound veloc- 
ity diminishes as the superfluid goes from the BCS to 
the unitary regime on increasing the strength of the at- 
tractions between the atoms, and that in both limits the 
speed of sound differs by a factor •v/3/5 from that in the 
superfluid three-dimensional gas. The numerical study 
of the superfluid gas dynamics in a weak axial confine- 
ment shows that the value of the speed of sound around 
in the central region of the trap is close to that in the 
fully quasi-onedimensional gas. However, in situations 
of experimental relevance, a density perturbation of or- 



4 



1 

T X T 
„ BB » J v 


\ T 




1 

- normal I— B—\ ' 
superfluid — G — l , 




i 


l/\/5 




i i i i 



1 10 10 2 10 3 10' 1 10 6 10 6 10 7 

|o|/oo 

FIG. 4: Sound velocity in an elongated trap as a function 
of \a\ (in units of the Bohr radius ao and in a log scale). 
Circles and squares are for a superfluid and for a normal Fermi 
gas, respectively. The error bars for the normal-gas results 
are from the standard deviation in the velocity distributions. 
The dot-dashed lines show the zero-sound, the BCS, and the 
unitary limit, from top to bottom. 

der 20% could give rise to upward shifts of the speed of 
sound by 10-12%. These effects could be observed in ex- 
perimental set-ups of trapped superfluid Fermi gases and 
thereby provide further insight in their superfluid state. 
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